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The mass spectra and electromagnetic decay rates of charmonium, bottomonium 
and B c mesons are comprehensively investigated in the relativistic quark model. The 
presence of only heavy quarks allows the expansion in powers of their velocities. All 
relativistic corrections of order v 2 /c 2 , including retardation effects and one- loop ra- 
diative corrections, are systematically taken into account in the computations of the 
mass spectra. The obtained wave functions are used for the calculation of radiative 
magnetic dipole (Ml) and electric dipole (El) transitions. It is found that relativistic 
effects play a substantial role. Their account and the proper choice of the Lorentz 
structure of the quark-antiquark interaction in a meson is crucial for bringing theo- 
retical predictions in accord with experimental data. A detailed comparison of the 
calculated decay rates and branching fractions with available experimental data for 
radiative decays of charmonium and bottomonium is presented. The possibilities 
to observe the currently missing spin-singlet S and P states as well as D states in 
bottomonium are discussed. The results for B c masses and decays are compared 
with other quark model predictions. 
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I. INTRODUCTION 



The investigation of the properties of mesons composed of a heavy quark and antiquark 
(bb, cc, cb) gives very important insight into heavy quark dynamics. Heavy quarkonia have 
a rich spectroscopy with many narrow states lying under the threshold of open flavour 
production. Excited states experience different decays among which there are radiative 
transitions to lower levels. The theoretical analysis shows that many properties of heavy 
quarkonia, including mass spectra and radiative decay rates, are significantly influenced 
by relativistic effects. Thus their inclusion is necessary for the correct description of the 
spectroscopy and the determination of quarkonium wave functions. Radiative decays are 
the most sensitive to relativistic effects. Some of these decays, which are forbidden in the 
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exact nonrelativistic limit (so-called hindered transitions) due to the orthogonality of initial 
and final meson wave functions, have decay rates of the same order as the allowed ones. 
In the relativistic description of mesons an important role is played by properties of the 
confining quark- ant iquark interaction, in particular its Lorentz structure. Thus comparison 
of theoretical predictions with experimental data can provide valuable information on the 
form of the confining potential. Such information is of great practical interest, since at 
present it is not possible to obtain the QQ potential in the whole range of distances from 
the basic principles of quantum chromodynamics (QCD). As it is well known, the growing of 
the strong coupling constant with distance makes perturbation theory inapplicable at large 
distances (in the infrared region). In this region it is necessary to account for nonperturbative 
effects connected with the complicated structure of the QCD vacuum. All this leads to a 
theoretical uncertainty in the QQ potential at large and intermediate distances. It is just in 
this region of large and intermediate distances that most of the basic meson characteristics 
are formed. 

At present, a vast set of experimental data is available on the masses and different decays 
of heavy quarkonia. However, not all states predicted by theory have been observed yet, 
while the others need confirmation. Such missing or unconfirmed states are present both in 
charmonium (2 1 S , , l 1 / 3 !, l 1 ^, 1 ^2,3) and bottomonium (spin-singlet and l P\ states, 
D states). The different possibilities for their experimental observation are proposed and 
widely discussed in the literature [1-6]. Radiative transitions from the spin-triplet levels with 
J PC = 1 to these states as well as their subsequent radiative decays play an important role 
in these proposals. The missing charmonium states can also be searched in B meson decays 
and identified by their radiative transitions [1, 2]. For this purpose, reliable relativistic 
predictions for the masses of these states and for the rates of radiative transitions involving 
them are necessary. 

The properties of the B c meson are of special interest, since it is the only heavy meson 
consisting of two heavy quarks with different flavour. This difference of quark flavours 
forbids annihilation into gluons. As a result, the excited B c meson states lying below the 
BD production threshold undergo pionic or radiative transitions to the ground pseudoscalar 
state which then decays weakly. There should be a rather rich set of such narrow states 
which are considerably more stable than corresponding charmonium or bottomonium states. 
The CDF collaboration [7] reported the discovery of the B c ground state in pp collisions at 
Fermilab. More experimental data are expected to come in the near future from new hadronic 
colliders. 

The purpose of this paper is to give a detailed analysis of mass spectra and radiative 
transitions in charmonium, bottomonium and B c mesons with the comprehensive account 
of the relativistic effects. This will allow one to get valuable information about the Lorentz 
structure of confining quark interactions from the comparison of obtained predictions with 
available experimental data. On the other hand, it will indicate the processes in which the 
missing states can be searched for. 

The paper is organized as follows. In Sec. II we describe our relativistic quark model. 
The expression for the heavy quark- ant iquark quasipotential with the account of relativistic 
(including retardation effects) and one loop radiative corrections is given in Sec. III. There it 
is applied to the calculation of the charmonium, bottomonium and B c meson mass spectra. 
In Sec. IV pseudoscalar and vector decay constants of the B c meson are calculated with the 
account of relativistic corrections and compared with other theoretical predictions. In Sec. V 
the relativistic expressions for the radiative transition matrix elements in the quasipotential 
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approach are given. They are used for the calculation of the decay rates of radiative Ml 
and El transitions in Sees. VI and VII, respectively. The role of relativistic effects in 
these transitions is investigated. Special attention is payed to the influence of the Lorentz 
structure of the quark potential on the relativistic corrections to decay rates. Pure vector 
and scalar potentials as well as their mixture are considered. The obtained results are 
compared with available experimental data, and the possibilities for searching the missing 
states in bottomonium are discussed. Finally, our conclusions are given in Sec. VIII. 



II. RELATIVISTIC QUARK MODEL 

In the quasipotential approach a meson is described by the wave function of the bound 
quark- ant iquark state, which satisfies the quasipotential equation [8] of the Schrodinger type 
[9] 

where the relativistic reduced mass is 

E l E 2 M 4 - (m\ - ml) 2 
m ~ E[+E~ 2 ~ 4AP ' (2) 

and Ei, E 2 are given by 

= M 2 -m 2 + m\ = M 2 - m? + m 2 

1 2M ' 2 2M U 

Here M = Ei + E 2 is the meson mass, mi^ are the quark masses, and p is their relative 
momentum. In the center of mass system the relative momentum squared on mass shell 
reads 

2 = [M 2 - (mi + m 2 ) 2 ][M 2 - (m t - m 2 ) 2 ] 
1 ] AM 2 ' 1 ] 

The kernel V(p, q; M) in Eq. (1) is the quasipotential operator of the quark- ant iquark in- 
teraction. It is constructed with the help of the off-mass-shell scattering amplitude, projected 
onto the positive energy states. Constructing the quasipotential of the quark-antiquark in- 
teraction, we have assumed that the effective interaction is the sum of the usual one-gluon 
exchange term with the mixture of long-range vector and scalar linear confining potentials, 
where the vector confining potential contains the Pauli interaction. The quasipotential is 
then defined by [10] 

V(p, q; M) = Ui(p)u 2 (-p)V(p, q; M)ui(q)u 2 (-q), (5) 

with 

V(p, q; M) = ^ SJ CV(k)7i72 + (k)rf r 2;/1 + vi nf (k), 

where as is the QCD coupling constant, is the gluon propagator in the Coulomb gauge 
D°°(k) = -^, D«(k) = -f 2 U* - ^) , £>« = J D» = 0, (6) 
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and k = p — q; 7^ and u{p) are the Dirac matrices and spinors 

u\p) 



e(p) + m 



1 
crp 



\ HP) \ e (p)+m 
with e(p) = \Jp 2 + m 2 . The effective long-range vector vertex is given by 



(7) 



IK 



r„(k) = 7^ + 2^^ fcV > 



(8) 



where k is the Pauli interaction constant characterizing the anomalous chromomagnetic 
moment of quarks. Vector and scalar confining potentials in the nonrelativistic limit reduce 
to 



V v {r) = (1 - e)Ar + B, 

V s (r) = eAr, (9) 

reproducing 

V COQ{ (r) = Vs(r) + V v (r) = Ar + B, (10) 

where e is the mixing coefficient. 

The expression for the quasipotential of the heavy quarkonia, expanded in v 2 /c 2 without 
and with retardation corrections to the confining potential, can be found in Refs. [10] and 
[11], respectively. The structure of the spin-dependent interaction is in agreement with 
the parameterization of Eichten and Feinberg [12]. The quasipotential for the heavy quark 
interaction with light antiquark without employing the expansion in inverse powers of the 
light quark mass is given in Ref. [13]. All the parameters of our model such as quark masses, 
parameters of the linear confining potential A and B, mixing coefficient e and anomalous 
chromomagnetic quark moment k are fixed from the analysis of heavy quarkonium masses 
[11] (see Sec. Ill) and radiative decays [14] (see Sees. V-VII). The quark masses m& = 4.88 
GeV, m c = 1.55 GeV and the parameters of the linear potential A = 0.18 GeV 2 and 
B = —0.16 GeV have usual values of quark models. The value of the mixing coefficient of 
vector and scalar confining potentials e — — 1 has been determined from the consideration 
of the heavy quark expansion for the semileptonic B — > D decays [15] and charmonium 
radiative decays [14]. Finally, the universal Pauli interaction constant k — —1 has been 
fixed from the analysis of the fine splitting of heavy quarkonia 3 Pj- states [10]. Note that 
the long-range magnetic contribution to the potential in our model is proportional to (1 + k) 
and thus vanishes for the chosen value of k — — 1. In the present paper we will take into 
consideration the retardation corrections as well as one-loop radiative corrections. 



III. HEAVY QUARKONIUM AND B c MESON MASS SPECTRA 

The heavy quark-antiquark potential with the account of retardation effects and one 
loop radiative corrections can be presented in the form of a sum of spin-independent and 
spin-dependent parts. For the spin-independent part we have [11] 

V SI ( r ) = JM5 + Ar + 5 _!M^1M 
3 r 3 2-7T r 
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+ 



+ 



2m\m 2 y \ 

4^«M 
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1 - e 



4 ay 
3 r 



P H y 



w 



2%) M 



m(/jr) 



1 



2m\m 2 4 y mf 



m 2 



P 2 - 
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r r J 
(P • r) 2 



- - 

4 Ira? 



+ 
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where 
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Bp 2 , 



14 2/ 7T 



(11) 



(12) 



«i = y - y n /> 



A) 



11 "a"'" 



Here n/ is the number of flavours and /x is a renormalization scale. Note that for a quantity 
quadratic in the momenta we use the Weyl prescription [16]: 

{f(rW}w = \{{f(r),P*}^}. 

For the dependence of the QCD coupling constant a s (fi 2 ) on the renormalization point 
/i 2 we use the leading order result 



47T 



A,ln(/i 2 /A 2 ) 



(13) 



In our numerical calculations we set the renormalization scale \x = 2m\m 2 j {m\ + ra 2 ) and 
A = 0.168 GeV, which gives a s = 0.314 for rai = ra 2 = m c (charmonium); a s = 0.223 for 
m i = m 2 = m b (bottomonium); and a s = 0.265 for rai = m c , m 2 = m b (B c meson). 

The spin-dependent part of the quark-antiquark potential for unequal quark masses with 
the inclusion of radiative corrections [17, 18] can be presented in our model as follows: 



V SB (r) =a L-S + b 



-(S 1 -r)(S 2 -r)-(S 1 -S 2 ) 



+ cSi-S 2 + dL.(Si-S 2 ), (14) 




3 12 lE 



y -3J + yln(/Jr) 



-3 ln( v /m 1 m 2 r) 



+ 



m\m<i 3 r 3 



1 + 



a s (/x 2 



7T 



1 

6 



12 
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+ie (y ~ ^ + y ln ^ r ) - \ HV m i m 2r) 




b = 



mi 

1 (4a s (/x 2 ) a s {fi 2 ) 



, rn 2 1/1 1 \ 2 , w . A 



m 2 / 



(15) 



3mim 2 ( 

—3 ln( v /mim 2 r 

4 J87ra s (/i 2 
3mim 2 \ 3 



7T 



)+(l + «) 2 (l-£)^}, 
a s (/i 2 



1 + 
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(16) 
<5 3 (r) 



a s (/i 2 



7T 



V 2 An(/xr) + 7 g \ [ ln (V m i m 2 r ) + 7g 
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+ (1 + k) 2 (1 -e)— , 
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(17) 



1 + 



a s (/i 2 



7T 
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where L is the orbital momentum and Si j2 , S = Si + S 2 are the spin momenta. For the equal 
mass case (mi = m 2 = m) the second order in a s contribution of the annihilation diagrams 



5c 



8^(/i 2 
3m 2 



1 -ln2)£ 3 (r) 



(19) 



must be added to the spin-spin interaction coefficient c in Eq. (17). 

The correct description of the fine structure of the heavy quarkonium mass spectrum 
requires the vanishing of the vector confinement contribution. This can be achieved by 
setting 1 + k — 0, i.e., the total long-range quark chromomagnetic moment equals zero, 
which is in accord with the flux tube [19] and minimal area [16, 20] models. One can see 
from Eq. (14) that for the spin-dependent part of the potential this conjecture is equivalent 
to the assumption about the scalar structure of confining interaction [21]. 

To calculate the heavy meson mass spectra with the account of all relativistic corrections 
(including retardation effects) of order v 2 /c 2 and one-loop radiative corrections we substitute 
the quasipotential which is a sum of the spin-independent (11) and spin-dependent (14) parts 
into the quasipotential equation (1). Then we multiply the resulting expression from the 
left by the quasipotential wave function of a bound state and integrate with respect to the 
relative momentum. Taking into account the accuracy of the calculations, we can use for 
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the resulting matrix elements the wave functions of Eq. (1) with the static potential 1 

V NR (r) = ~^± + Ar + B. (20) 
3 r 

As a result we obtain the mass formula 



b 2 (M) 
where 



W + (a) (L • S) + (b) ( A( Sl • r) (S 2 • r) - (S 1 • S 2 ) ) + (c) (S x • S 2 > + d(L ■ (S 1 - S 2 )> , 

(21) 



W = (Vsi) + (22) 

The first term on the right-hand side of the mass formula (21) contains all spin-independent 
contributions, the second and the last terms describe the spin-orbit interaction, the third 
term is responsible for the tensor interaction, while the forth term gives the spin-spin inter- 
action. The last term is not zero only for the unequal mass case mi ^ m 2 and leads to the 
mixing of triplet and singlet meson states with the total angular momentum J equal to the 
orbital momentum L. 

In Table I the calculated charmonium mass spectrum is compared with experimental 
data. For meson states we use the notation n 2S+l Lj, where n = n r + 1 and n r is the radial 
quantum number. Our predictions agree with PDG [22] data within a few MeV. Our model 
correctly reproduces both the position of the levels and their fine and hyperfine splittings. In 
this Table we also give the recent Belle Collaboration data [23] on pseudoscalar rj c (l ^So) and 
rj' c (2 ^o) states observed in B decays. The measured mass of the ground spin-singlet state 
7] C {1 1 Sq) is in good agreement with world averages and predictions of our model, while the 
radially excited r]' c {2 l S$) state lies considerably higher than previous experimental indications 
and most of the theoretical predictions. If these data are confirmed, it will be difficult to 
accommodate such a small hyperfine splitting ss 32 MeV (almost four times smaller than 
117 MeV splitting for the ground state) in the framework of the quark model. 2 

Our prediction for the mass of h c (l : Pi) is consistent with the data from the Fermilab 
Experiment E760 [25] on pp — * /i c (3526) — > tt°J/^ which, however, need confirmation. 
The same is true for the indication of a 1 3 D 2 state with mass 3836 ± 13 MeV in 7r ± N — > 
J/^vr+TT^+anything [26]. 

For the calculation of the bottomonium mass spectrum it is also necessary to take into 
account additional one-loop corrections due to the finite mass of the charm quark [27-30]. 
We considered these corrections within our model in Ref. [31] and found that they give 
contributions of a few MeV and are weakly dependent on the quantum numbers of the 
bottomonium states. The one-loop correction to the static QQ potential in QCD due to the 
finite c quark mass is given by [29, 31] 

AV(r, m c ) = —-^ [ln(V^m c r) + lE + E^^m.r)} , (23) 
y nr 



1 This static potential includes also some radiative corrections. The remaining radiative correction term 
with logarithm in Eq. (11), also not vanishing in the static limit, is treated perturbatively. 

2 The position of the ^'(2 3 Si) can in principle be influenced by the nearby threshold of the open charm 
production. 
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TABLE I: Charmonium mass spectrum (in GeV). 



State 


Particle 


1 neory 


Experiment 


2S+lr 
n °^ Lj 






t "> t "\ r\ Tool t ) n Tool 

FDij [22\ Belle [23J 


i 1 in 

1 S 


Vc 


2.979 


2.9797 2.979 


i So 
l a Si 


J/* 


3.096 


3.09687 


l 6 P 


XcO 


3.424 


3.4151 


-I ^7-» 
l 6 Pl 


Xcl 


3.510 


3.51051 


1 


Xc2 


3.556 


3.55618 


-i 1 r> 


h c 


3.526 




o lc 

2 


Vc 


o COO 


3.594 3.654 






3.686 


3.68596 


i 3n 




3. r9o 


3. r o9y 


1 3 D 2 




3.813 




1 £>3 




3.815 




1 




3.811 




2 3 P 


XcO 


3.854 




2 3 Pi 


Xcl 


3.929 




2 3 P 2 


Xc2 


3.972 




2^1 




3.945 








3.991 




3 3 5i 




4.088 


4.040 6 


"This value from Ref. [24] is included only in 


the PDG Listings. 





^Mixture of S and D states. 



where 

7# = 0.5772 is the Euler constant and a = 5.2. The resulting bottomonium mass spectrum 
with the account of this correction is given in Table II. We found that the small shift of the 
QCD parameter A from our previous [11] value 0.178 to 0.168 GeV (with all other parameters 
remaining fixed) allows us to get a good fit to the bottomonium mass spectrum with the 
account of finite charm quark mass corrections. The difference between the theoretical and 
experimental data is less than 3 MeV. Very recently CLEO Collaboration presented [32] 
the first evidence for the production of the triplet Y(1P) state in the four photon cascades 
starting from T(3S). In Table II we give their preliminary result for the mass of T(1I?2) 
state which is consistent with our prediction. 

For the mesons consisting of quarks with different flavours (mi ^ m 2 ), such as the B c 
meson, the coefficient d (18) in the spin-dependent part of the quark potential (14) is not 
equal to zero. This results in the mixing of singlet and triplet P states with J = 1, 

nPl' = n x Pi cos 9 n p + n 3 Pi sin 9 n p, 

nPl = -n x Pi sin 6> nP + n 3 Pi cos G nP , (24) 
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TABLE II: Bottomonium mass spectrum (in GeV). 



btate 


Particle 


1 neory 




Experiment 








i ~) i "\ r \ Tool 

FUG [22\ 


r \ r / \ Tool 

CLEO [32J 


1 


Vb 


9.400 








T 


9.460 


9.46030 




1 3 P 


Xbo 


9.863 


9.8599 






Xbi 


9.892 


9.8927 




1^2 


Xb2 


9.913 


9.9126 




1 Pi 


7 

h b 


9.901 






n In 


i 

% 


9.993 






2 cS 5i 


T' 


10.023 


10.02326 




1 3 Di 




10.153 






1 3 £>2 




10.158 




10.162 


1 3 7^1 




10.16z 






1^2 




10.158 






2 3 P 


Xm 


10.234 


10.2321 




2 3 Pi 


X&i 


10.255 


10.2552 


10.2556 


2 3 P 2 


Xfe2 


10.268 


10.2685 


10.2688 


2^1 




10.261 






3^0 




10.328 






3 3 5i 




10.355 


10.3552 





and of D states with J = 2, 

n.02' = n ^ cos + n 3 D 2 sin #„d, 

n£>2 = -n^ sin 6 nD + n 3 D 2 cos 9 nD . (25) 
For the B c meson the values of the mixing angles in our model are 

e lP = 0.357, e 2P = 0.405, e 1D = -0.627. (26) 

In Table III we compare our model predictions for the mass spectrum of the B c meson 
with other quark model results [33-36] . We see that the differences between the predictions 
in most cases do not exceed 30 MeV. The only exceptions are masses of ID states, which 
are 50 — 70 MeV heavier in our model. The fine and hyperfine splittings are also consistent 
with each other. All these predictions for the ground state pseudoscalar B c and vector B* 
meson masses satisfy the bounds found by Kwong and Rosner [37]: 

6.194 GeV < M Bc < 6.292 GeV 

and 

6.284 GeV < M B , < 6.357 GeV. 

c 

In Ref. [38] the ground state B c mass was evaluated in perturbative QCD. Experimental 
data [22] at present are available only for the B c ground state and have large error bars 
M Bc = 6.4 ± 0.4 GeV. 
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TABLE III: B c meson mass spectrum (in GeV). 



State 
n 2S+1 Lj 


Our 


[33] 


[34] 


[35] 


[36] 


l^o 


6.270 


6.264 


6.253 


6.286 


> 6.2196 


l 3 ^i 


6.332 


6.337 


6.317 


6.341 


> 6.2786 


1 3 P 


6.699 


6.700 


6.683 


6.701 


> 6.6386 


1P1 


6.734 


6.730 


6.717 


6.737 


> 6.7012 


lPl' 


6.749 


6.736 


6.729 


6.760 


> 6.7012 


1 3 P 2 


6.762 


6.747 


6.743 


6.772 


> 6.7347 


2^0 


6.835 


6.856 


6.867 


6.882 




2 3 5i 


6.881 


6.899 


6.902 


6.914 




l 3 P>i 


7.072 


7.012 


7.008 


7.019 




ID2 


7.077 


7.009 


7.001 


7.028 




1D2' 


7.079 


7.012 


7.016 


7.028 




1 3 D 3 


7.081 


7.005 


7.007 


7.032 




2 3 P 


7.091 


7.108 


7.088 






2P1 


7.126 


7.135 


7.113 






2Pl' 


7.145 


7.142 


7.124 






2 3 P 2 


7.156 


7.153 


7.134 






3^0 


7.193 


7.244 








3 3 5i 


7.235 


7.280 









In the following sections we apply the masses and wave functions of T and B c mesons 
for the calculation of their decay constants and decay rates. 

IV. PSEUDOSCALAR AND VECTOR DECAY CONSTANTS OF THE B c 

MESON 

The B c meson and its first excitations which lie below the BD threshold are stable against 
strong decays, since they cannot annihilate into gluons. They can decay via electromagnetic 
and pionic transitions into the lightest pseudoscalar ground state B c . The significant contri- 
bution to the B c total decay rate comes from the annihilation of the c quark and b antiquark 
into the vector boson W + which decays into a lepton and a neutrino or a quark- antiquark 
pair. The weak annihilation decay rate is determined by the pseudoscalar constant of the 
B c meson. 

The decay constants fp and fy of the pseudoscalar (P) and vector (V) mesons param- 
eterize the matrix elements of the weak current between the corresponding meson and the 
vacuum. In the case of the B c meson they are defined by 

(0|6 7 At 7 5 c|P(K)) = ij P K\ (27) 
(0\h^c\V(K,e)) = f v M v e^ (28) 

where K is the meson momentum, P corresponds to the pseudoscalar B c and V to the vector 
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B* mesons, and My are the polarisation vector and mass of the vector B* meson. 

In the relativistic quark model the decay constants can be expressed through the meson 
wave function $py(p) in the momentum space and are given by 



fp,\ 




e b (p) +m b \ 
2e 6 (p) ) 



1/2 



[e c (p) +m c ][e b (p) + m b ] 



Spyip), (29) 



with Ap = — 1 and Ay = 1/3. In the nonrelativistic limit p 2 /m 2 — > these expressions for 
decay constants give the well-known formula 



/^ = /^ = J^l*/>.r(0)|. (30) 



py 

where ^py(0) is the meson wave function at the origin r = 0. 

The calculated values of the pseudoscalar and vector decay constants of the B c meson in 
our model using the relativistic formula (29) are displayed in Table IV. They are compared 
with the ones calculated using the nonrelativistic expression (30) and other predictions of 
the nonrelativistic quark models [33, 35], QCD sum rules [34] and lattice NRQCD [40]. 
We see that inclusion of relativistic corrections reduces the pseudoscalar decay constant fp c 
by 20% and produces the difference between vector and pseudoscalar decay constants of 
approximately 70 MeV. The calculated values of these decay constants are consistent with 
lattice [40] and QCD sum rule [34] predictions. 



TABLE IV: Pseudoscalar and vector decay constants (fp = /b c , fv = /bj) of the B c meson (in 
MeV). 



Constant 


rel 


NR 


[33] 


[34] 


[35] 


[40] 


!b c 


433 


562 


500 


460 db 60 


517 


420 ± 13 


Ib* 


503 


562 


500 


460 ± 60 


517 





V. RADIATIVE TRANSITIONS IN HEAVY QUARKONIA AND B c MESONS 

To determine the rates of radiative decays (B — > A + 7) it is necessary to calculate the 
matrix element of the electromagnetic current between the initial (B) and final (A) meson 
states. In the quasipotential approach such matrix element has the form [41] 

04 = / ^^^ P(P) ^ (P ' <*>*B Q (q), (31) 

where T^(p,q) is the two-particle vertex function and ^a,b are the meson wave functions 
projected onto the positive energy states of quarks and boosted to the moving reference 
frame. The contributions to Y come from Figs. 1 and 2. The contribution r^ 2 ^ is the 
consequence of the projection onto the positive-energy states. Note that the form of the 
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FIG. 1: Lowest order vertex function I^ 1 ) corresponding to Eq. (33). Radiation only from one 
quark is shown. 
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FIG. 2: Vertex function corresponding to Eq. (34). Dashed lines represent the interaction 
operator V in Eq. (5). Bold lines denote the negative-energy part of the quark propagator. As on 
Fig. 1, radiation only from one quark is shown. 

relativistic corrections resulting from the vertex function explicitly depends on the 
Lorentz structure of the QQ-interaction. Thus the vertex function is given by 



where 
and 



r M (p,q) = r( 1 )(p,q) + r( 2 )(p,q) + ---, 

r«(p,q) = e 1 u 1 (p 1 ) lfl u 1 (q 1 )(2n) 3 5(p 2 - q 2 ) + (1 <- 2), 



It^q) = e 1 u 1 (p 1 )u 2 (p 2 ){v(p 2 - q 2 ) ^ ^ 



ei(fci) + ei(gi) 



+7im ^l! (fcl l , 7i°V( P2 - q 2 )}«i(?i)« 2 (ft) + (1 - 2). 



Here ei j2 are the quark charges, k x = p x — A; = q x + A; A = P — Q; 

a^Hp) = e(p) " (m l°A 7 ° (7 ' p)) , C (p) = 4v^\ 



Hp) 



(32) 
(33) 



(34) 



and 



Pi,2 = e 1)2 (p)|^±5> w (^y, 
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9i,2 = e 1)2 (g)^±£n«( PB )g\ 

where are three four-vectors given by 

" ( ' , "W = {f^« + m|tm)}' b =v^^- « = 1 - 2 - 3 - 

= (Eg, Q) and = (-Ea, P) are four-momenta of initial and final mesons. 

It is important to note that the wave functions entering the current matrix element (31) 
cannot be both in the rest frame. In the initial B meson rest frame, the final A meson is 
moving with the recoil momentum A. The wave function of the moving A meson f aa is 
connected with the wave function in the rest frame ^! aq = ^ a by the transformation [41] 

*aa(p) = D\ /2 (R^)D\ l2 (R^ AO ( V ), (35) 

where R w is the Wigner rotation, L A is the Lorentz boost from the rest frame to a moving 
one, and the rotation matrix D l l 2 (R) in the spinor representation is given by 

(J ^D[^RYj = S-\ Pl , 2 )S(A)S(p), (36) 

where 

/ e(p) + m ( a p \ 

5(p) = V^^( 1 + ^)T^J 

is the usual Lorentz transformation matrix of the four-spinor. 

To calculate the radiative transition matrix element we adopt the following procedure. 
We substitute the vertex functions and given by Eqs. (33) and (34) in the decay 
matrix element (31) and take into account the wave function transformation (35). The 
resulting structure of this matrix element is rather complicated, because it is necessary to 
integrate both over d 3 p and d 3 q. The 5 function in expression (33) permits us to perform one 
of these integrations and thus this contribution can be easily calculated. The calculation of 
the vertex function r^ 2 ^ contribution is more difficult. Here, instead of a 5 function, we have 
a complicated structure, containing the QQ interaction potential in the meson. However, 
we can expand this contribution in powers of the heavy quark velocities v 2 /c 2 and then use 
the quasipotential equation in order to perform one of the integrations in the current matrix 
element. It is easy to see that the vertex function contributes already at the first order 
of the v 2 /c 2 expansion. 

We consider two main types of radiative transitions: 

a) Magnetic dipole (Ml) transitions which go with the spin flip of the quark (AS = 1, 
AL = 0) and thus the initial and final states belong to the same orbital excitation but have 
different spins. Examples of such transitions are vector to pseudoscalar (n 3 Si — > n' + 7, 
n > n') and pseudoscalar to vector (n x Sq — > n' 3 S± + 7, n > n') meson decays. 

b) Electric dipole (El) transitions in which the orbital quantum number is changed (AL = 1, 
AS = 0) and thus the initial and final states belong to different orbital excitations but have 
the same spin. Examples of such transitions are n 3 Si — > n' 3 Pj7 (n > n') and n 3 Pj — > n' 3 S'i7 
(n > n') decays. 
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VI. RADIATIVE Ml TRANSITIONS 



A. Ml decay rates 



The radiative Ml transition rate is given by [14] 



r(B 



A + -f) = —{2J 1 + 1)\M B a\ , where u 
Sit 



M 2 B -M\ 
2Mr ' 



(37) 



Mb and Ma are the initial and final meson masses, J' is the total angular momentum of 
the final meson. The matrix element of the magnetic moment is defined by 



M BA = -- 



d 

dA 



x (A |J(0)| B) 



A = P - Q, 



(38) 



A=0 



where (A |J M (0)| B) is the matrix element of the electromagnetic current between initial (B) 
and final (A) meson states with momenta Q and P, respectively. 

After inserting the vertex functions and from Eqs. (33) and (34) in the decay 
matrix element (31) with the account of the wave function transformation (35), we carry 
out the expansion in inverse powers of the heavy meson masses Mr^, which are large due to 
the presence of two heavy quarks Mb, a ~ + m Q'- Then we calculate the matrix element 
of the magnetic moment operator (38) and get 
(a) for the vector potential 



M 



v 



d 3 p 



ei 



(l-e)(l + 2/c)[px [ffixp] 



(2tt)3 AV ^2 Cl (p) 
'1 -e)(l + K)\p x [cr 2 x p] 



2e 1 (p)[e 1 (p) + mi] 



ei(p)[e2(p) + m 2 ] 



+ 



C2(P) 

' M B 
1 



l + (l-e) 



M B - e 1 (p) - e 2 (p)' 



d_ 

^ dp 



P x 



2Mb 

(b) for the scalar potential 



P x 



cr 2 



d 3 P t i \ 
*a(p) 



ei 



, e i(p) + m i ti{v) + m i , 



e^p) + e 2 {p) - M B 



>*b(p) + (1<-2), (39) 



(2k) 



2ei(p) 



x 



e 2(p) 
<Ti 



2Mp 



M f 



P x 



P x 




P x 



01 



e l(p) 

e [p x [o-i x p]] 

2ei(p)[ei(p) + mi 



ei(p)+mi e 2 (p)+m 2 



>* B (p) + (1^2). (40) 



Note that the last terms in Eqs. (39), (40) result from the wave function transformation (35) 
from the moving reference frame to the rest one. It is easy to see that in the limit p/m — > 
the usual nonrelativistic expression for the magnetic moment follows. 

Since we are interested in radiative transitions between S state (vector and pseudoscalar) 
mesons, it is possible to evaluate spin matrix elements using the relation (ai) = — (cr 2 ). 
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Then, taking into account that both quarks are heavy (Q and Q'), we further expand 
Eqs. (39), (40) in inverse powers of the heavy quark mass itlq up to the second order 
corrections to the leading contribution and get 

(a) for the purely vector potential (e = 0) 

2m Q { 3m z Q 3 m Q \m Q m QI J 6M B \m Q m Q > J J 

(41) 

(b) for the purely scalar potential (e — 1) 

Ms = j3_L_ ^ + M(i _ + -L)U Q ~cn 

2m Q { m Q 2m Q \mQ> 3m Q J 6M B \ m Q m Q' J J 

Here (• • •) denotes the matrix element between radial meson wave functions. For these 
matrix element calculations we use the meson wave functions obtained calculating their 
mass spectra. 

For the quarks of the same flavour (rag = rag/, eg = —eg/) and k — — 1 these expressions 
reduce further [14] 

(a) for the purely vector potential (e = 0) 

Mr^JSLk-W^LX (43) 

2m Q [ 3m 2 Q 3M B m Q J 

(b) for the purely scalar potential (e — 1) 

^ S 2 2mg{ 3 m Q + 3m Q 3M B mg}' 



B. Results and discussion 



The resulting Ml radiative decay rates of charmonium, bottomonium and B c are pre- 
sented in Tables V-VII. In these Tables we give predictions both for allowed (n 3 Si — > 
n' + 7, n' — n) and hindered (n > n') decays. For the calculation of allowed decay rates 
we use expanded expressions (41)-(44). For the hindered transitions, which are strongly 
suppressed in the nonrelativistic limit due to orthogonality of the initial and final state wave 
functions, relativistic effects are decisive. Thus for their calculations we use unexpanded 
expressions (39) and (40). In Tables V-VII we present the photon energy u, the decay rates 
calculated discarding all relativistic corrections T NR , as well as using relativistic expressions 
for purely vector T v , for purely scalar V s and for the mixture (9) of vector and scalar poten- 
tials T with e — — 1. Note that in all these calculations we use the relativistic wave functions 
found calculating the meson mass spectra in Sec. III. 

The Ml radiative decay rates are very sensitive to relativistic effects. Even for allowed 
transitions relativistic and nonrelativistic results differ significantly. An important example 
is the decay J/\I/ — > ^7. It is well known that the nonrelativistic predictions for its rate are 
more than two times larger than the experimental data. As we see from Table V, the inclu- 
sion of the relativistic corrections for purely scalar or purely vector potentials do not bring 
theoretical results in agreement with experiment. For the purely scalar potential the decay 
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TABLE V: Radiative Ml decay rates of charmonium. For decays involving rj' c we give in parenthesis 
the results obtained using the recent value [23] of its mass. The values T exp are taken from Ref. [22]. 



Decay 


to 


pNR 


T v 


T s 


L 


pcxp 




MeV 


keV 


keV 


keV 


keV 


keV 


J/* r/ c7 


115 


2.73 


1.95 


3.13 


1.05 


1.13 ±0.35 


*' Vol 


91(32) 


1.26(0.055) 


0.85(0.037) 


0.71(0.031) 


0.99(0.043) 


0.6 - 3.9 a 




639 


0.23 


0.61 


0.35 


0.95 


0.84 ±0.19 


rf c -► J/*7 


463(514) 


0.26(0.36) 


0.70(0.95) 


0.37(0.51) 


1.12(1.53) 




°This value from Ref. [24] needs confirmation and is included only in the PDG Listings. 






TABLE VI: Radiative Ml decay rates of bottomonium. 




Decay 




pNR 


T v 


T s 


r 


B 




MeV 


eV 


eV 


eV 


eV 


(io- 4 ) 


T -> ?7 b 7 


60 


9.7 


8.7 


12.2 


5.8 


1.1 


T' - ^7 


33 


1.6 


1.45 


1.50 


1.40 


0.32 


T" - <7 


27 


0.9 


0.8 


0.8 


0.8 


0.30 


T' -> 7? 6 7 


604 


1.3 


3.4 


1.3 


6.4 


1.5 


< ^7 


516 


2.4 


6.3 


2.5 


11.8 




T" r? b7 


911 


2.5 


6.2 


3.1 


10.5 


4.0 


< - T 7 


831 


5.8 


14.3 


7.1 


24.0 




T" - ^7 


359 


0.2 


0.6 


0.1 


1.5 


0.57 


^ - T' 7 


301 


0.4 


1.1 


0.2 


2.8 





rate even increases by 15%. On the other hand for the purely vector potential relativistic 
effects decrease the decay rate by 25%, but such decrease is not enough: the theoretical 
result still deviates from experimental data by more than 2a. 3 Only for the mixture of 
vector and scalar potentials (9) we get the necessary decrease of the decay rate which brings 
theory in agreement with experimental data for the J/^f — > r\ c ^ decay rate. For the hindered 



TABLE VII: Radiative Ml decay rates of the B c meson. 



Transition 




-pNR 


T v 


V s 


r 


L [33] 


r [34] 


r [35] 




MeV 


eV 


eV 


eV 


eV 


eV 


eV 


eV 


l 3 ^ -l 1 ^ 


62 


73 


48 


66 


33 


135 


60 


59 


2 3 Si-2 1 S 7 


46 


30 


24 


32 


17 


29 


10 


12 


2 3 S 1 -l 1 ^ 


584 


141 


412 


398 


428 


123 


98 


122 


2 1 So^2 3 S 1 7 


484 


160 


471 


454 


488 


93 


96 


139 



3 This is compatible with the estimate that relativistic effects can give contributions of order of 20-30% in 
charmonium. 
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decay — > r? c 7 the decay rate calculated for the mixture of vector and scalar potentials is 
also in good agreement with experiment while the rates for pure potentials (especially the 
scalar one) are lower than the experimental value. In Table V we give predictions for decays 
involving the first radial excitation of the pseudoscalar state rj' c (2 ^o) as well. Since there 
are two contradicting experimental measurements of its mass we calculated the rates using 
both values. The results obtained using recent Belle value [23] are given in parenthesis. 

In Table VI predictions for Ml decay rates and branching fractions of bottomonium are 
given. Since the hyperfine splitting in bottomonium is predicted to be small (around 60 
MeV, see Table II) the photon energies and hence decay rates of allowed Ml transitions are 
very small. This is one of the main reasons why no spin-singlet S-wave levels rjbin^o) have 
been observed yet. Our results show that relativistic effects for the favored mixture of vector 
and scalar potentials further decrease the allowed Ml decay rates. 

Recently it was argued by Godfrey and Rosner [4] that hindered transitions could be 
more favorable for discovering the rjb- Their analysis of different quark model predictions 
showed that most nonrelativistic models favor the rj^lS) production from T(2S) decays, 
while the account of relativistic corrections makes prospects for discovering rjb(lS) in T(3S) 
radiative decays comparable to those in T(2S*) decays. Our present relativistic consideration 
of these decays supports this observation. Indeed we see from Table VI that relativistic 
effects significantly increase rates of hindered transitions. The hindered decay T"(3S) — > 
7]b(lS) , ~f has the largest branching fraction 4.0 x 10 -4 , which is almost 2.7 times larger than 
the T'(2S') — > 77;,(l t S')'7 decay branching fraction. Very recently CLEO Collaboration [42] 
searched for r] b (lS) in such hindered Ml transition from T"(3S). No evidence of T"(3S) — > 
r/fe(15')7 transitions was found and rather strict upper limits on the branching fraction were 
set: B(T"(3S) — > rjb(lS)^) < 6 x 10~ 4 , which rule out many previous phenomenological 
predictions reviewed in Ref. [4]. Our model result for the branching fraction of this decay is 
below but rather close to this experimental upper limit. 

In Table VII we give predictions for decay rates of Ml radiative transitions of the B c 
meson in our model in comparison with previous nonrelativistic quark model analysis [33- 
35]. We see that relativistic effects play an important role in B c meson Ml radiative decays. 
They reduce the rates of allowed decays and increase the rates of hindered transitions. The 
largest rates are predicted for the latter decays which are increased by relativistic effects 
almost by the factor of 3 and thus they are an order of magnitude larger than the rates of 
allowed Ml transitions. 



VII. RADIATIVE El TRANSITIONS 



A. El decay rates 

The radiative El transition rate is given by [43] 

T(B^A + 1 ) = ^-\B BA \ 2 , where u = M \~f A , (45) 

Mb and Ma are the initial and final meson masses. The matrix element of the electric dipole 
moment Dba is defined by 



B BA = -i-^(A\J {0)\B) 



A = P - Q, (46) 

A=0 
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where (A 1^(0)1 B) is the matrix element of the electromagnetic current between initial (B) 
and final (A) meson states with momenta Q and P, respectively. 

We substitute expressions (31)-(34) in the definition of the electric dipole moment (46) 
and take into account the relativistic transformation of the wave function (35). Then, 
discarding some terms of order v A /c A and higher, we get the following expressions for the 
electric dipole moment Dba [43] (indices 1,2 are changed to Q,Q') 

(a) for the purely vector potential 

n - ( d3p * M P U €Q,{p) 9 x p] (-, £q(p) 

+ 2[ Mb -^-^] \_ [^xp] 1 ( p )-(Q„<y), (47) 
m Q J 2M B [e Q/ {p) +m Q ,\ J 

(b) for the purely scalar potential 



(2ir)' "^'^ M B dp 2i Q (p){e Q (p) + m Q ]\ M B 
_ 2[M fl - tg (ri- fQ ,(ril \ _ [ CTg ,xpl | M _, (4g) 
m Q / 2M B [e Q ,(p) +m Q ,\ J 

The operator id /dp in Eqs. (47), (48) corresponds in the coordinate space to the operator 
r. All other terms in these equations are relativistic corrections. Thus in the nonrelativistic 
limit the standard expression for the electric dipole moment is recovered. 

It is easy to see that there are three different structures with respect to the orbital 
variables in Eqs. (47), (48): r, [(ctq + ctq^xp] and [(<tq — <tq/) xp]. Thus the matrix element 
of the electric dipole moment for the electromagnetic transition nJMLS — > n'J'M'L'S' + 7 
can be presented in the form 

B vs = (n'J'M'L'S'\A(p 2 )r - B v > s (p 2 )[S x p] - C v ' s (p 2 )[(S Q - S Q ,) x p]\nJMLS), (49) 

where functions A(p 2 ), B v ' s (p 2 ) and C v ' s (p 2 ) up to order p 2 jm 2 Q^Q^ are given by 

,/ 2n _ e Q m Qf -e Q >m Q e Q m Q - e Q >m Q > p 2 
{P ] ~ M B + M B 2m Q m QI ' (5Uj 

B v (p 2 ) = J9- (1 + 2(M ~ mQ ~ mQ,) 7p2 p2 ^ 



4m 2 Q \ m Q 4m 2 , m 2 QI ) 

e Q' ( 1+ 2(M-m Q -m Q >) 7p 2 p 2 



4m 2 Q , \ m Q/ 4m 2 Q , m\ 

e Q + e Q > \ 1 1 p 2 ( 1 1 



AM 



B 



m Q m Q , 4 \m 3 Q m\, J _ 
2{M-m Q -m Q ,) 7p 2 p 2 



(51) 



c v (p 2 ) = -pv 1 + 

4m 2 Q \ m Q 4m 2 Q m 2 



+-V 1 + 



Q'J 

2(M -rriQ-mey) 7p 2 p 2 



4m 2 ,, V m Q , Am 2 QI m% 
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eg + e Q > 



4M f 



1 



1 



ttlq rriQi 



P_ 

4 



1 1 
— H - 

Q 



m 



(52) 



B s (p 2 " 



4m 2 Q 



' x _ 2(M -m Q - rngQ | p 2 | 



m Q 



4m 2 



q rn Q m Q > 



eg' 



4m 2 Q , . 
eg + e Q / 



i _ 2(M - m Q - mgQ p 



+ 



4M f 



1 



1 



4m 2 Q , m Q mQi 



mQ vtiqi 



4 



4m^ 



i _ 2(M - m Q - m Q /) p 



2 



(53) 



+ 



4m 2 Q m Q m QI 



4m Q , 



^ 2{M-m Q -m Q ,) | p 



+ 



™g' 



4m 2 Q , m Q m QI 



eg + e Q / 
4M R 



1 1 

+ 



™g rn QI 



4 



1 1 



(54) 



The last structure in Eq. (49) proportional to [(Sq — Sq/) xp] leads to the spin-flip transitions. 
It vanishes for the cc and bb mesons, consisting of the quark and antiquark of the same flavour 
(itlq = rnqr, e Q = —e Q >) since in that case C v ' s (p 2 ) = 0. The functions A(p 2 ), B v ' s (p 2 ) then 
simplify and coincide with the ones found previously [43] 



A(p 2 ) = 
B V ( P 2 ) = 
B S (P 2 ) = 



+ 



eg 

eg 
2m 2 



' 2m Q 
M B ^ M B m Q J ' 

2{M B -2m Q ) llp r 



1 + 



4m3 



' x _ 2{M B - 2m Q ) + 5p 2 



4mj 



In the case of the B c meson C v,s {p 2 ) ^ and thus relativistic corrections lead to spin-flip 
transitions (S' = S ± 1) but only for decays involving mixed states nPl, nPV (24) or 
nD2, nD2' (25). For all other transitions the spin- flip correction vanishes due to momenta 
relations (see Eq. (57) below). 

Using the Wigner-Eckart theorem and relations for matrix elements of the tensor operator 
between coupled functions, one can rewrite Eq. (49) in the form 



sJ'+J+L+S-M' 



J' 1 J 

—M' i M 



(2j + i)(2j' + 1; 



X 





J' 




[{'; 


L 


n 



- X (J', L', J, L)8 SS ' ± i (n'L'\\C(p 2 )p\ \nL 



(55) 
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where 



r)(J', L', J, L, S) = (-1) L+J+1 V6[S(S+1)(2S+ I)} 1 ' 2 



X (J',L',J,L) = (-1) 



L+L' 



L' J' S 
J L 1 

/ 2 / J L 1 
2 J' + 1 i 1 1 V 



-i \ J S L 

J' S L' }, (56) 

1 1 1 

(57) 



Here 



J' 1 J 
-M' g M 



L' J' S 
J L 1 



and < 



J S L 
J' S V 
1 1 1 



are 3j-, 6j- and 9j-symbols, (n'L'H • • • \ \nL) are reduced matrix elements. 

The total El decay rate of the nJLS state is obtained by summing the decay rates (45) 
over all possible values of M' for a fixed value of M. The resulting expression is 



F v ' s (nJLS -> n'J'L'S' + 7) = -a to 3 1 5 SS >C 1/2 (J\ L', J, L, S) 



(L' - L) lR n , L , A(p 2 )r 



R 



nL 



_r,(J',L',J,L,S) / n > L >\ lB v, Sip 2 MnL 
Jmax(L', L) 



S S s'±iX(J',L',J,L)(n'L'\\C v ' s (p 2 )p\\nL) , (58) 



where 



C(f, L', J, L, S) = max(L', L)(2J' + 1) 



V J' S 
J L 1 



The reduced matrix elements can be expressed through the usual matrix elements over radial 
wave functions i? n ^(r): 

(a) for the transitions between P and S states 



(n'0|H|nl) = -(R' n , s \R nP ), 

[n'0\\ V 2 p\\nl) = (RZ s \R nP )+2^K, s 
(b) for the transitions between D and P states 



RnP ) — 2 ( R' , 



R 



nP 



(59) 



(n'l|H|n2) = -V2[(R n ,p 



RnD ) — (R n 'p\RnD) , 



n'l||p 2 p||n2) = -V2[{RZ P \R nD ) + {K, 



-6 ( R' n , P 



1 

r 

RnD ) + 6 ( R n 'P 



R 



nD 



R 



nD 



(60) 



where the prime means differentiation of R with respect to r. 
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TABLE VIII: Radiative El transition rates of charmonium. For decays involving rf c we give in 
parenthesis the results obtained using the recent value [23] of its mass. 



ci y 


1 1 I 


-pNR 


r v 

L 




p 


pexp fool 




MeV 

1VJ.C V 


keV 


keV 


keV 


keV 


keV 


z, D\ > i -T07 


9^Q 


M 7 


04.U 


44 n 

44. U 


9fi 3 


o« i _|_ q o 

Z*U . X ZIZ O . 


o3c, . I 3p 


1 71 


4:4. cf 


3D 1 


38 3 
oo.o 


99 Q 


o _i_ q pi 

—■•J.Z, ZIZ O.U 


9 3c v, — > 1 3 p n0/ 


1 9R 


^0 Q 
ou.y 


99 Q 


98 1 


1 S 9 


on 4 4- 9 ^ 

ZU.4 ZTZ Zr.O 


9 1< ?n _> 1 ip.'V 

' "JO > 1 r Y) 


fi^n 9&^ 


O.U I (J / J 


9f41 ^ 


fi 9(41 "1 


9(41 1 




1 3 Pn — > 1 3< ? 1 'V 
1 -TQ ► 1 Ol7 


OUO 


1 fil 
±u± 


1 ^1 


1 S4 


1 91 

X At X 


IUiJ ZIZ ou 


1 3 P, — > 1 3< ? 1 'V 

i .ri ► i 017 




ooo 


98^ 


OUO 




iijl ZIZ O-L 


1 3 Pn — > 1 ^i-V 

i x^2 ' 1 '-'17 




4zLR 


QflQ 


9Q9 


397 


OOcJ ZIZ tJZi 


1 ri — > 1 Oo7 




79^ 


ODU 


ODU 


ODU 




l 3 ^i 1 3 P 7 


361 


423 


344 


334 


355 




1 3 Z>i -► 1 3 Pi7 


277 


142 


127 


120 


135 




1 3 Di 1 3 P 27 


234 


5.8 


6.2 


5.6 


6.9 




1 3 D 2 -► 1 3 P l7 


291 


297 


215 


215 


215 




1 3 D 2 -► 1 3 P 2 7 


248 


62 


55 


51 


59 




1 3 £> 3 ^1 3 P 27 


250 


252 


163 


170 


156 




1 ^2 1 l P Xl 


275 


335 


245 


245 


245 





B. Results and discussion 

The results of numerical calculations of charmonium El radiative decay rates using 
Eqs. (56)-(60) are presented in Table VIII. For calculations of photon energies u we used the 
experimentally measured masses of charmonium S and P states. 4 For masses of D states 
we used our model predictions from Table I. We give predictions for decay rates calculated 
in the nonrelativistic limit T NR , for relativistic decay rates with pure vector T v and scalar 
V s potentials as well as for the mixture (9) of vector and scalar potentials T with e — — 1. 
As in the case of Ml decay rates calculations, we use the relativistic wave functions in our 
numerical analysis. 

The results presented in Table VIII show that relativistic effects play an important role 
in El decays of charmonium. The most sensitive to the relativistic corrections are decays 
ty'(2S) — > XcjO-P) + 7- Their account leads to the considerable reduction of the decay 
rates. The rates for the vector potential are reduced more significantly than for the scalar 
one. As a result, there arises an approximately twofold reduction of decay rates for the 
mixture of vector and scalar potentials with the value of mixing parameter e — — 1, bringing 
theoretical predictions in good agreement with experimental data. The large influence of 
relativistic corrections originates from the fact that the zero of the 25* wave function is close 
to the maximum of the IP wave function. This results in a reduction of the leading order 
decay matrix element (lP|r|25'|). Therefore relatively small relativistic corrections produce 



4 For decays involving ij' c (2S), as in the case of Ml transitions, we use both experimental values of its mass, 
giving a prediction for the recent Belle value [23] in parenthesis. 
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such a large effect. This observation is confirmed by the calculations of the Xcj(l-P) ~~ * 
J/ty + 7 decay rates. Here both initial and final states do not have zeros and relativistic 
contributions have usual values and lead to an approximately 25% reduction of the decay 
rate. All theoretical predictions are in nice agreement with data. In Table VIII we also give 
predictions for El decay rates of charmonium ID states. At present only the 1 3 D\ state 
is experimentally observed. This state is considerably broader, since it lies above the DD 
threshold. The observed state ^(3770) can also have a significant 2 3 S± state admixture [44]. 
If we consider it to be a pure D state, then using its measured total decay rate, we get the 
following predictions for the El radiative decay branching fractions: 

B{1 3 £>i -> 1 3 P 7) ~ 1-5%; B{\ 3 D X -> 1 3 P l7 ) « 0.6%; B(l 3 P>i -> 1 3 P 27 ) « 0.03%. 

On the other hand, the tensor 1 3 D2 and 1 1 D 2 charmonium states are under the threshold 
of open charm production, since their masses are slightly below the DD* threshold, and the 
decay of these states to DD is forbidden by parity and angular momentum conservation. 
Thus El radiative transitions are the main decay channels of these states. 

The calculated decay rates of El radiative transitions in bottomonium are presented 
in Tables IX and X. The influence of relativistic effects in bottomonium is considerably 
less than in charmonium. The contribution of relativistic corrections does not exceed 10% 
almost for all decays. The only exceptions are decays T"(3S) — > XbjO-P) + 7> where the 
leading contribution is substantially reduced due to the significantly different number of 
zeros in initial 3S and final IP wave functions. For all S — > P + 7 transitions we find good 
agreement of our model predictions with experimental data. 

The comparison of the theoretical predictions for the radiative decays of P states of 
bottomonium XbijiP) with the experimental data is complicated by the fact that the total 
decay rates of these states are not measured yet. Experiment gives only branching fractions 
B = T[xb(nP) — > T(n'S) + 7] /r total (n > n') for these decays. Thus for this comparison 
it is necessary to get theoretical predictions for the total decay rates of Xb{ n P)- The main 
decay channels of the bottomonium P states are inclusive strong decays to gluon and quark 
states and radiative decays. The strong decays were extensively studied in the literature 
[45, 46] including leading-order QCD corrections. The two-gluon annihilation rates of 3 Po 
states with the account of relativistic corrections are given in our model [47] by 



r( 3 P - 99) 



5<X 



3M 2 



+ 1 



d 3 p m 2 Q 
(2tt) 3 Mp 

—)( 

2e Q (p)J \ 



P 



In 



eg(p) +P 



c q(p) ln e Q(p)+P X 



2 - 



P ' ^q{p)-P / 



>(P) 



(61) 



To include leading-order QCD corrections [46] this expression should be multiplied by (1 + 
10.0o; s /7r) for Xbo an d (1 + 10.2a s /7r) for x'w The corresponding expression for two-gluon 
annihilation rates of 3 P 2 states reads as [47] 



n 3 p 2 



99) = 



_8o|_f r d 3 p m Q e Q (p) 
5M 2 \ J (2tt) 3 Mp 



2 + 



P 



£q(p) 



ln £q(p) +P _ £q{p) ln eg(p) + P' 



x 1 + 



m Q 



£q(p) 
2p 2 



P 



2eg ip) [eg (p) + m Q ] J 3e Q (p) [e Q (p) + m Q ] 



p ^q(p)-p, 



<P) 
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TABLE IX: Radiative El transition rates of bottomonium. 



Decay 




pNR 


T v 


V s 


r 


r exp [22] 




MeV 


keV 


keV 


keV 


keV 


keV 


2 3 5i -► l 3 P o7 


162 


1.65 


1.64 


1.66 


1.62 


1.67 ±0.37 


2 3 5i -► l 3 Pi7 


130 


2.57 


2.48 


2.51 


2.45 


2.9 ±0.6 


2 3 5i -► 1 3 P 27 


109 


2.53 


2.49 


2.52 


2.46 


3.08 ±0.6 


2 1 S -l 1 Pi7 


98 


3.25 


3.09 


3.09 


3.09 




1 3 P ^ 1 3 ^17 


391 


29.5 


30.6 


31.3 


29.9 




1 3 Pi -► 1 3 ,Si7 


422 


37.1 


37.0 


37.4 


36.6 




1 3 P 2 1 3 ^17 


442 


42.7 


39.8 


39.3 


40.2 




1 X P! - 1 ^07 


480 


54.4 


52.6 


52.6 


52.6 




3 3 5i ^2 3 P 07 


123 


1.65 


1.51 


1.52 


1.49 


1.42 ±0.25 


3 3 ^i -► 2 3 Pi 7 


100 


2.65 


2.43 


2.45 


2.41 


2.97 ±0.43 


3 3 Si ^2 3 P 27 


86 


2.89 


2.69 


2.71 


2.67 


3.0 ± 0.45 


3 1 5o->2 1 Pi 7 


73 


3.07 


2.78 


2.78 


2.78 




3 3 5i 1 3 P 07 


484 


0.124 


0.040 


0.054 


0.027 




3 3 Si ->■ l 3 Pi7 


453 


0.307 


0.097 


0.134 


0.067 


0.041 ±0.029 a 


3 3 5i 1 3 P 27 


433 


0.445 


0.141 


0.195 


0.097 


0.064 ±0.045 a 


3 ^o^l^n 


427 


0.770 


0.348 


0.348 


0.348 




2 3 P ^2 3 S l7 


207 


11.7 


11.1 


11.2 


11.0 




2 3 P X -► 2 3 5 l7 


230 


15.9 


14.8 


14.8 


14.7 




2 3 P 2 -► 2 3 S"i 7 


243 


18.8 


16.7 


16.6 


16.7 




2 ip x _> 2 1 5 o7 


262 


23.6 


21.4 


21.4 


21.4 




2 3 P ^ l 3 S l7 


743 


7.36 


7.58 


8.41 


6.79 




2 3 Pl -► l 3 ,S l7 


764 


8.01 


7.90 


8.33 


7.49 




2 3 P 2 l 3 S l7 


776 


8.41 


7.61 


7.20 


8.02 




2^1 ^1^07 


820 


9.9 


9.36 


9.36 


9.36 





"Needs confirmation. 




ln £ ^ )+P + 3 

CQ{P)-P V P 

eg(p) ln e Q(P) +P ] + 



p2 



P £q(p)-PJ m Q [e Q {p) + m Q ] 




34 (P) 

p2 



H 



'(P) 



M \ 
m Q J 

(62) 



with leading-order QCD corrections [46] accounted by the factors (1 — 0.1a s /7r) for Xb2 and 
(1 + 1.0a s /7r) for x'b2- Here <f>p(p) is the P state radial wave function in momentum space. 
In the nonrelativistic limit p/mg — > and M — > 2mg, Eqs. (61) and (62) reduce to the 
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TABLE X: Radiative El transition rates of bottomonium involving D states. 







pNR 


T v 


r 5 


r 




MeV 


keV 


keV 


keV 


keV 


1 3 Pi — ► 1 3 Pn^ 


280 


24.2 


23.4 


23.4 


23.4 




256 


12.9 


12 7 


12 7 


12 7 


1 3 Di — > 1 ^o'v 


235 


67 


69 


69 


69 


1 3 Do — > 1 3 Pi^V 


969 


94 8 


93 3 


93 3 


93 3 


1 3 Do — > 1 3 Ptt 


241 


6.45 


6.35 


6.35 


6.35 


1 3 P/q — ► 1 3 Po^ 


244 


26 7 


24.6 


24.6 


24.6 


1 1 Do — > 1 X Pi^ 
J- x J 1 / 


254 


SO 2 


28.4 


28.4 


28.4 


9 3 Pn — > i 3 ni^ 
z j^o ' 1 - L/ i y 


O-L 


1 1 7 


1 16 


1 16 


1 1 7 

J. . J. 1 


2 3 Pi -► l 3 L>i7 


104 


0.62 


0.61 


0.60 


0.615 


2 3 Pl ^ 


98 


1.56 


1.56 


1.56 


1.56 


2 3 P 2 _> l 3 D l7 


117 


0.036 


0.034 


0.034 


0.035 


2 3 P 2 _> l 3 D 27 


111 


0.454 


0.446 


0.443 


0.449 


2 3 P 2 l 3 D 37 


108 


2.34 


2.36 


2.37 


2.35 


2 ^ -► 1 1 Z? 27 


102 


2.62 


2.43 


2.43 


2.43 



known expressions [46]: 



rft. M )^, (63) 



r( 3 P 2 - a,) = 8a -'ft (0)l • (64) 



The calculation of radiative and relativistic corrections to the annihilation decay rates of 
3 Pi and l P\ states is a very complicated problem which has not been solved yet. Thus to 
estimate their hadronic decay rates we use the tree-level nonrelativistic expressions [46] 

r( 3 Pi -?? + </) = ^l\K P m 2 Hm Q (r)), (65) 
20a? _ |2 . 



T( 1 P 1 -> ^) = ^— A- |K P (0)| 2 ln(m Q (r», (66) 



T( 1 P 1 - ^ + 7 ) = f eJ-r( 1 P 1 -> «?«?«?). (67) 

5 q a s 

From the relativistic consideration of the decays of 3 P 0j 2 states, which shows that relativistic 
effects give ~ 10% contributions to the bb decay rates, we can expect that these formulae 
give a reasonable estimate of the corresponding decay rates. For numerical calculations of 
hadronic decay rates of Xb states we use a s = 0.18 obtained from the experimental ratio of 
T(T — > ggj)/T{T — > ggg) [5]. The calculated partial decay rates and branching fractions 
for IP and 2P states of the bottomonium are compared with available experimental data 
in Table XI. There we give both PDG [22] averages and very recent CLEO [32] data. We 
see that our predictions for the branching fractions for radiative transitions of Xb( n P) to 
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T(n'S) states are in good agreement with experiment. The only discrepancy (of 2.5er) is 
the CLEO value for B(xbi(2P) ~^ T(2S) + 7) transition which is approximately two times 
larger than our model prediction and PDG value. The CLEO Collaboration [32] measured 
also two photon decays of T(3S) via Xb(^-Pj) states. They report the branching fraction for 
T(3S) — > XbO-Pj) + 7 T(IS') + 77 transitions summed over all the J states: 

B(T(3S) -> T(1S) + 77) = (2.14 ± 0.22 ± 0.21) x 10 -3 . (68) 

Using our model results in Tables IX-XI for corresponding decay rates, we get 

B(T(3S) - X w(lP) + 7 ^ T(15) + 77) = 1-4 x lO" 5 , 
B(T(3S) - XbiiXP) + 7 ^ T(IS-) + 77) = 9.97 x 10" 4 , 
£(T(3S) -> Xi 2 (2P) + 7 -> T(15) + 77) = 9.96 x 10' 4 , 

and the sum over all the J states is equal to 

B(T(3S) -> T(1S) +77) = 2.04 x 10~ 3 

in accord with CLEO data (68). 

The CLEO Collaboration [32] presented recently the first evidence for the production 
of the triplet T(1D) states in the four photon transitions 3S — > 2P + 7 — > IP + 77 — > 
IP + 777 — > I5 + 7777 — > e + e _ + 7777. The measured product branching fraction for these 
five decays is equal to (3.3±0.6±0.5) x 10 -5 . In Table XII we give the theoretical predictions 
for the branching fractions of such four photon decays in our model and in the recent 
quark model analysis of Godfrey and Rosner [6]. In general, both theoretical predictions 
are consistent. In the analysis of Ref. [6], the dominant decay is 3 3 Si — > 2 3 P 1 + 7 — > 
1 3 D 2 + 77 — > 1 3 -Pi + 777 — > 1 3 S± + 7777 — > e + e~ + 7777, while in our model the above and 
3 — > 2 3 P 2 + 7 — > 1 ^ + 77 — > 1 3 P2 + 777 — >■ 1 + 7777 — >■ e + e _ + 7777 transitions 
have almost the same rate and dominate. In the last line we give the sum of all these decay 
channels. Both theoretical predictions agree with CLEO measurement. 

An important problem of quarkonium physics is the search for bottomonium spin singlet 
states rjb(n ^o) and hb(n x Pi). In the previous section we discussed the possibilities of finding 
rji, in radiative Ml decays. From Table XI we see that these two states can be discovered 
simultaneously. Indeed the radiative decay hj,(l 4 Pi) — > ^(1 ^o) +7 with the photon energy 
of 480 MeV is the main decay channel of (the branching fraction of this decay exceeds 
50%). Thus production of a few h b states, e.g., through T"(3S) -> h b (l ^tttt or T"(3S) -> 
h b (l 1 Pi)7r decays, which branching fractions are predicted to be about 0.1-1% [48], will give 
a good possibility of finding r] b . 

In Tables XIII and XIV we compare our results for the El radiative decay rates of 
the B c meson with other quark model predictions [33-35]. Comparison of the calculations 
using relativistic T and nonrelativistic T NR formulae for decay rates shows that relativistic 
corrections do not exceed 20% in B c meson El decays. Most of the theoretical predictions 
for El transitions between P and S states of B c mesons given in Table XIII are compatible 
with each other. The largest differences occur for decays involving PI and PI' states which 
are the mixtures of spin singlet and spin triplet states (24) due to different mixing angles 
used by the authors. Note that for such transitions there are additional relativistic spin-flip 
contributions to decay rates (58) proportional to x{J' ■> L' , J, L) (57) which are specific only 
for B c mesons. In general, our predictions are closer to the ones of Ref. [34]. 

In Table XIV we present the El radiative decay rates of B c mesons where either the 
initial or final state is a D wave state. Here we find rather large variations in theoretical 
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TABLE XI: Partial decay rates and branching fractions for IP and 2P states of bottomonium. 



Level 


Decay 


r (keV) 


B (%) 


PDG [22] 


CLEO [32] 


1 3 P 


55 


653 


95.6 








l 3 Si + 7 


29.9 


4.4 


< 6 




l 3 Pi 


qq + g 


57 


60.9 








l 3 Si + 7 


36.6 


39.1 


35 ±8 




1 3 P 2 


99 


109 


73 








l 3 Si + 7 


40.2 


27 


22 ±4 




l^i 


999 


36 


40.1 








99 + 1 


1.2 


1.3 








1 % + 7 


52.6 


58.6 






2 3 P 




431 


95.8 








2^ + 7 


11.0 


2.4 


4.6 ±2.1 


2.59 ±0.92 ±0.51 




l 3 5i + 7 


6.8 


1.5 


0.9 ±0.6 


< 1.44 




1 3 #i + 7 


1.17 


0.3 






2 3 Pi 


39 + 5 


50 


67.2 








2 3 S 1 + 7 


14.7 


19.8 


21 ±4 


41.5 ±1.2 ±5.9 




l 3 5i + 7 


7.5 


10.0 


8.5 ± 1.3 


11.6 ±0.4 ±0.9 




l 3 #i + 7 


0.6 


0.8 








1 3 D 2 + 7 


1.6 


2.2 






2 3 P 2 




76 


73.4 








2^ + 7 


16.7 


16.1 


16.2 ±2.4 


19.3 ±1.1 ±3.1 




l 3 5i + 7 


8.0 


7.7 


7.1 ± 1.0 


7.0 ±0.4 ±0.8 




l 3 #i + 7 


0.04 


0.03 








1 3 £> 2 + 7 


0.5 


0.5 








l 3 P/ 3 + 7 


2.4 


2.3 






2 x Pi 


555 


31.5 


47.9 








55 + 7 


1.1 


1.7 








2% + 7 


21.4 


32.5 








1 % + 7 


9.4 


14.3 








1 X D 2 + 7 


2.4 


3.6 







predictions. The main reason of these distinctions is the difference in values of D state 
masses, which for some states reaches 70 MeV (see Table III). Since 2P and ID states of 
the B c are rather close, such a difference significantly influences the energy of the emitted 
photon and thus the decay rates. For decays involving the mixed spin singlet and spin triplet 
states PI, PI' (24) and D2, D2' (25) the additional relativistic spin-flip contributions (58) 
are important, especially for transitions where both initial and final states are mixed states. 
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TABLE XII: Predicted branching fractions of four photon decays of T(3S) involving ID states 
corresponding to 3S — ► 2P + 7 — ► ID + 77 — ► IP + 777 — ► IS 1 + 7777 — ► e + e~ + 7777 transition. 
The last line gives the sum of all decay channels. 



2 3 Pj state 


1 3 D j state 


1 3 Pj state 


gour 

(xl0~ 6 ) 


B [61 
(xl0~ 6 ) 


2 3 P 2 




l 3 Po 


15.2 


7.8 






1 3 P> 


0.7 


0.3 






1 3 P 


3.9 


2.7 




l 3 Dl 


1 3 P 9 


0.0 


0.0 






l 3 Pi 


0.1 


0.1 






l 3 Pn 


0.0 


0.0 




1 3n 


1 3p 


z.y 


z.o 






l 3 Pl 


15.5 


20.1 




l 3 £>i 


1 3 P 2 


0.1 


0.1 






l 3 Pl 


2.4 


3.3 






1 3 P 


0.5 


0.4 


2 3 P 


l 3 ^! 


1 3 P 2 


0.0 


0.0 






l 3 Pl 


0.5 


0.3 






1 3 P 


0.0 


0.0 


all 


all 


all 


41.8 


37.6 



VIII. CONCLUSIONS 

In this paper the mass spectra and radiative Ml and El decay rates of charmonium, 
bottomonium and B c mesons were calculated in the relativistic quark model based on the 
quasipotential approach in quantum field theory. Special attention was devoted to the role 
of the relativistic effects in these processes. Since both quarks in the considered mesons 
are heavy, the v/c expansion was applied. In the mass spectra calculations retardation as 
well as one-loop radiative corrections were taken into account. We also included the one- 
loop correction due to the finite charm quark mass to the bottomonium mass spectrum. It 
was found that this correction is rather small [31] and its inclusion allows one to obtain an 
even better fit of the bottomonium excited states with the slightly shifted value of QCD 
parameter A. The calculated mass spectra of charmonium and bottomonium agree with 
the experimental data within a few MeV. Comparison of our results for the B c meson mass 
spectrum with previous calculations showed that different predictions for ground state and 
low excitations agree within 30 MeV. 

The pseudoscalar and vector decay constants of B c meson were calculated using the 
relativistic wave functions obtained during the mass spectrum calculations. It was found that 
relativistic effects reduce these constants by approximately 20% and produce the splitting 
between them of about 70 MeV. 

It was shown that relativistic effects play a significant role in radiative decays of mesons. 
Their form strongly depends on the Lorentz structure of the quark-antiquark interaction. 
The most sensitive are radiative Ml decays, where even for allowed transitions they sig- 
nificantly influence predictions for the rates. An important example is the J/ty — > rjcj 
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TABLE XIII: Radiative El transition rates of the B c meson. 



T~)poa,v 




pNR 


r v 


V s 


r 


r [331 


r [341 


r [351 




MeV 


keV 


keV 


keV 


keV 


keV 


keV 


keV 


1 3p _> j 3c ^ 


355 


75.5 


74.1 


81.3 


67.2 


79.2 


65.3 


74.2 


1P1 1 35 l7 


389 


87.1 


82.9 


87.1 


78.9 


99.5 


77.8 


75.8 


1P1' — »■ 1 3 S*i7 


405 


13.7 


12.6 


11.6 


13.6 


0.1 


8.1 


26.2 


1 3 P> — ► 1 ■y 


416 


122 


105 


102 


107 


112.6 


102.9 


126 


1P1 — > 1 ^rTV 


447 


18.4 


16.3 


14.4 


18.4 


0.0 


11.6 


32.5 


1P1' _». 1 lQ n ™ 


463 


147 


134 


136 


132 


56.4 


131.1 


128 




181 


5.53 


4.36 


5.00 


3.78 


7.8 


7.7 


9.6 


2 3 S*i — > IPI7 


146 


7.65 


5.98 


6.98 


5.05 


14.5 


12.8 


13.3 


2 3 Si — > IPI'7 


130 


0.74 


0.62 


0.61 


0.63 


0.0 


1.0 


2.5 


2 35 _> 1 3p > ^ 


118 


7.59 


5.99 


6.86 


5.18 


17.7 


14.8 


14.5 


9 — > IPI^ 


101 


1.05 


0.94 


0.86 


1.02 


0.0 


1.9 


6.4 


2 — > 1 Pl'<7 
w (J ~ ii x y 


84 


4 40 


3 77 

O.I 1 


3 82 


3 72 


5 2 


15 9 

X<J . U 


1 3 1 


2 3 P ^2 3 S l7 


207 


34.0 


28.8 


28.4 


29.2 


41.2 


25.5 




2P1 -► 2 3 5i7 


241 


45.3 


37.6 


37.3 


37.9 


54.3 


32.1 




2P1' -► 2 3 5i7 


259 


10.4 


8.45 


7.86 


9.07 


5.4 


5.9 




2 3 P 2 2 3 5 l7 


270 


75.3 


58.9 


60.6 


57.3 


73.8 


49.4 




2P1 -► 2 ^07 


285 


13.8 


11.1 


10.5 


11.7 




8.1 




2P1' -► 2 ^07 


303 


90.5 


73.2 


73.8 


72.5 




58.0 





transition, which is overestimated by a factor of more than two if the nonrelativistic approx- 
imation is used. It is argued that the inclusion of relativistic corrections for a pure scalar 
or vector confining potential is not enough to bring theoretical predictions in accord with 
experiment. Only for the specific mixture of these potentials (9) with the mixing coefficient 
e = — 1, the agreement can be obtained. For other decay rates this mixing of scalar and 
vector potentials also gives the best results. The hindered Ml transition rates are domi- 
nated by relativistic contributions and are significantly enhanced by them. The comparison 
of the allowed and hindered Ml rates in bottomonium shows that the latter provide better 
opportunity for the search of the missing pseudoscalar rjb state of the bottomonium. 

The analysis of radiative El transitions showed that the form of relativistic corrections 
is less dependent on the Lorentz structure of the quark potential than in the case of Ml 
transitions. However, for some decays, e.g., — > XcJl, the consideration of the mixed 
(9) vector and scalar potentials (with the same value of e — —1) is important for bringing 
decay rates in accord with experimental data. In general, all our predictions for radiative 
decay rates and branching fractions of charmonium and bottomonium agree with measured 
values. In the case of the B c meson radiative El decays an important additional relativistic 
correction to decay rates which causes the flip of the quark spin was found. This contribution 
(57) to the radiative El decay rate (58) is specific only for transitions involving mixed states 
nPl, nPV (24) or nD2, nD2' (25) of B c and is caused by the difference of the c and b quark 
masses. Finally, a comparison of various quark model predictions for the radiative Ml and 
El decay rates of B c has been performed. These radiative transitions along with pionic ones 
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TABLE XIV: Radiative El transition rates of the B c meson involving D states. 



Dprav 




pNR 


T v 


r s 


r 


r f33i 

1 [OOJ 


r [341 




MeV 


keV 


keV 


keV 


keV 


keV 


keV 


2 3 Po — ► 1 3 Di -7 


84 


0.035 


0.031 


0.027 


0.035 


0.2 


0.1 


2 3 Po — > 1 D2^ 


79 


285 


245 


0.222 


269 


3 2 

O . Zj 


1 5 

X . '.J 


2 3 Po — > 1 /52'^ 


77 


1 39 


0.114 


116 

\J • X X \J 


113 




5 


2 3 Po — > 1 3 n,^ 


75 


2 OS 


1 64 


1 69 
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are the main decay channels of the low lying excitations in the B c meson. 
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